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Intensity Problem in the Deflection of Cosmic Rays 
in the Solar Magnetic Field 


By Teiichiro YAGI 


Physical Institute, Nagoya University 


(Read May, 1951, Received Oct. 1, 1956) 
Abstract 


Assuming uniform parallel beam of cosmic rays at infinity is deflected in 
the general magnetic field of the sun, cosmic-ray intensities at each point of the 
earths’ orbit were calculated. Problem was treated only in the equatorial plane 
for simplicity. As the intensities change with season, the solar diurnal variation 
averaged over all season never tends to zero due to the heliomagnetic deflection 


of cosmic rays. 


We will assume that there exists a anisotropic distribution of cosmic rays in a 
space outside the solar system, due to any origin. This origin, for instance, may be 
the rotation of our galaxy, well known as Compton-Getting effect... Because of the 
deflection of cosmic-ray particles by the sun’s general magnetic field, apparent solar 
diurnal variation of cosmic rays will show an annual variation not only in its phase, 
but also its amplitude. Thus, even if this apparent solar diurnal variation will be 
averaged over one year, it will not vanish and remain as mean solar diurnal variation. 
In order to study this possibility, following calculation was carried out. Deflection of 
cosmic rays in the solar magnetic field was already calculated in our previous paper,” 
which will be denoted as [1] hereafter, and 
this work is continuation of [1]. 

Now, problem was treated only in the 
equatorial plane, under the assumption that 
the sun’s and the earth’s magnetic. dipole 
moment are parallel to each other. In Fig. 
1, M, is the sun’s dipole moment, which is 
directing towards the back of paper perpen- 
dicular to the plane of paper, that is, towards a 
the negative direction of Z-axis. Circle C is Fig. 1 Circle C det thet oer a! Genie and 
the earth’s orbit and P is the position of ee Mehee ay ana ee 
earth, which is designated by radius 7, and . 
angle y, counted positively clockwise from Y-axis. Now, we assume that there exists 
a uniform parallel beam of cosmic rays parallel to Y-axis at infinity. Furthermore, 
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we take y-axis perpendicular to Y-axis at sufficiently far distance from the earth and 
x-axis perpendicular to the direction of incidence of cosmic rays at P, and we consider 
that particles existing in the region ry~y+4y at infinity are incident into the region 
0~4x at P. Then, if we denote cosmic-ray intensities per unit time per unit area 
perpendicular to the direction of incidence of cosmic rays at infinity and at P as j(r) 


and j(¢,) respectively, then hare 
I(ys)dt dx dZ=](7)dt dr dZ. 


If we take the direction of x-axis in a definite way and j(y)=1, then 


¢ Or 
Pate ee ag. ; (1) 


F 0 0 
Hey=\ 22 ee a Or dr, 


dx | 0g, dx Or, dx 


where lengths are measured in St6rmer’s unit and y is two times 7,, which appears 
in [1]. 
Or 


0 ’ ee é 
and 2” were computed by differentiating the curve ¢,(7,, 7%:), which was 


0¢, Or, 
obtained at the study of [1] and is shown in Fig. 4a and 4b of appendix 1. So, of 
course, the accuracy of computation is rather low. Intensities 7(y,) thus obtained is 
shown in Fig. 2, where full line and dotted line are two branches with different 
directions of incidence of cosmic rays as already mentioned in [1]. From this, it is 
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Here, we will see how the direction of incidence 


of cosmic rays changes on x-axis. In Fig. 3, angle of 
incidence @, is counted positively eastwards from PO’. 
(In the paper [1], 8, was counted positively eastwards 
from PO’ for outer orbit and from. PO for inner 
orbit, respectively.) About the angle’ a from x-axis 
to the direction of incidence of cosmic rays, it is 


easily seen that 4a=40,—4¢,. Thus, 


da _d0, dy, _ dé, cosé, 


Se se eta: Ae = Fig. 3 Arc C is the earth’s 
a orbit and P is the position of 
da ~ 40,—cos0,° Ax (2) earth. Axes and symbols used 
Ys in the paper are shown. 
On the other hand, from sin 0, sigh wncicy 
Ae ee 
: d(sin 6, d : 
Tee fee eu ae E mt ot +- z sin 0, \« ae C3) )) 
dx dx i ee %, 


Thus, except the case that ue is considerably large, 4a is the order of oe If we 
take fifty times earth’s radius for 4x, which is considered to be region for cosmic rays 
being effectively influenced by the earth’s magnetic field, ae is 2.1x10-* and so Ja is 
negligible small. Thus, cosmic rays are incident into effective region of earth’s 
magnetic field in a parallel beam approximately. Furthermore, if we take usual 
dimension of cosmic-ray meter for 4x in the case of deflection of cosmic rays by the 
earth’s magnetic field, da is negligible small in this case also and cosmic rays are 
incident upon the erath in a parallel beam approximately. Thus, in the case that 
deflections of cosmic rays by the sun’s and the earth’s magnetic field are taken into 
account together, intensity distribution of cosmic rays on the earth in a certain season 
is such as the part corresponding to outer orbit in Fig. 2, but its relative magnitude 
changes with season in such a manner as shown in Fig. 2 also. 

Next, at the limit of increasing intensity in Fig. 2, j(g,). tends to « due to 


increasing of ae towards «. Although this divergence is merely mathematical and 


; eS tae, Aal dy ae oe aa 
integral intensity \ IPs) dz 2\i S| ge | dy is finite, the intensity is very strong at 
this point. But, at the same time, the angle 0, of incidence of cosmic rays changes 
very largely near this point, as known from equation 3. Intensity problem in such 
case cannot be treated simply and was omitted out of the present consideration. 
However, it is only small region of y, to be excluded for this reason. 

Now, the dependence of angle 0,(y,) of incidence of cosmic rays upon ¢, is given 
by Fig. A2 of appendix 1. From Fig. 2 and Fig. A2, we can express the intensity of 
cosmic rays at each point of the earth’s orbit with its direction of incidence expressed 
in solar time on harmonic dial. They are shown in Fig. 4a, 4b and ie There, two 
full lines correspond to two branches with different directions of incidence of cosmic 
rays as seen in Fig. 2 and dotted line is vector sum of these two branches. ~The 
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position on the earth’s orbit is shown by angle ¢, written on the curves. When we 
consider cosmic rays of vertical incidence on the earth, OY’-axis and OY’-axis ought 
to be taken instead of OY-axis in the cases of M,=10%4 and M,=2.5x10%3 gauss cm’, 


respectively. 


Fig. -4a For r,=4 Stormer of cosmic-ray 
energy, intensity of cosmic rays at each point 
of earth’s orbit with its direction of incidence 
are reproduced on harmonic dial. In the case 
that cosmic rays incident vertically on the 
earth are considered, axis OY’ and OY” are to 
be taken for o-hour axis instead of OY in 
the case of M,=10%4 and 2.5108 gauss cm‘, 
respectively. Two full lines correspond to two 
branches with different directions of incidence 
of cosmic rays as seen in Fig. 2 and dotted 
line is vector sum of these two branches. 
Vector A, B, and C are resultant vectors of 
two branches mentioned above and their vector 
sum. 
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| APPENDIX eee 
1. In Fig. Ala and Alb, angle ie t 
 ¢s is plotted against 71, taking” Ts 5 *r 
for parameter. The parts AeB-C a Tht 
of curves correspond to outer orbit. : Be 
‘In Fig. A2, angle 6, is plotted 4 ber 


taking r, for parameter. | Aa re. 
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Diurnal Variation of Vertical Cosmic Rays 
——Narrow Total and High Energy Components 


By Teiichfro YAGI and Hiroyuki UENO 


Physical Institute, Nagoya University 


(Read May 6, 1955 and May 14, 1956; Received Oct. 1, 1956) 
Abstract 


It is well known that anomalous solar diurnal variation Sp of cosmic rays 
appears at the time of geomagnetic storms. We confirmed the existence of Sp 
about the cosmic-ray intensities of total component and high energy component 
penetrating 90cm lead absorber measured with narrow angle telescope directed 
vertically. Amplitude of S, of the total component is the largest and the time of 
maximum intensities of that of the high energy component is the latest. The phases 
of Sp are different with five different kinds of cosmic-ray meters in our country. 
This can be understood qualitatively by the deflections of cosmic rays in the 
earth’s magnetic field. Year to year variation of phase of Spy and long-time 


variation of amplitude of mean diurnal variation were investigated. 
1. Introduction 


Since 1949, continuous observation of the intensities of cosmic ray total component 
has been carried out in our laboratory, using a narrow angle telescope directed 
vertically. Furthermore, since 1952, we have measured continuously the intensities of 
the high energy component penetrating 90cm lead absorber, with the same narrow 
angle telescope directed vertically. 

When geomagnetic storm occurs, the amplitude of diurnal variation of cosmic 
rays increases and the time of its maximum intensities advances. This phenomena 
were already found by Y. Sekido and S. Yoshida (1950) and D.W.N. Dolbear and H. 
Elliot (1951).2 We investigated in what manner this phenomena happen for the 
cosmic-ray intensities of high energy component penetrating 90cm lead absorber. In 
connection with this result, some features of anomalous diurnal variation of cosmic 


rays were reexaminend. 
2. Apparatus 


The apparatus is almost the same as described in our previous paper,” but some 
alterations are as follows. The telescope consists of 4 trays, having 6 Geiger-Miiller 
counters in each, as shown in Fig. 1. The effective length of individual counter is 
28cm and the effective diameter is 4,5cm. Between the bottom tray and the tray 
above it, 90cm lead absorber is inserted. The numbers of 3 fold coincidences (a, b, c) 
and (b, c, d) are recorded in every ten minutes, through scale-of-4 circuits respectively. 
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The maximum zenithal angles of incident cosmic-ray parti- 
cles causing these 3 fold coincidences are both 12 degrees, 
their averaged numbers per hour are 1744 and 87x4, 
and their statistical fluctuations of hourly values are 3.7% 
and 5.4%, respectively. The minimum momentum of 
cosmic-ray particles causing 3 fold coinciences (b, .c dias 
1,3 Gev/c for » mesons, owing to their penetration of 90cm 
lead. Also, the numbers of 2 fold coincidences (b, c) are | 
recorded in every minute, through a scale-of-16 circuit. 
The corresponding values of them are 40 degrees, 618 x 16 
particles per hour, and 0.98%, respectively. In the following, 
we will denote ZT (narrow angle, total), H (narrow angle, 
high energy) and W (wide angle, total) for 3 fold coinci- 
dences (a, b, c), (b, c, d) and 2 fold coincidence (b, c), 


respectively, in abbreviation. 
3. Barometer effect and mean diurnal variation 


3—1. Barometer effect The barometric pressure 
was observed at Nagoya District Meteorological Observa- 
< Piel pront view a tory, 16 times a day till December, 1952, and 8 times a day 
apparatus thereafter. The heights of isobar levels of upper atmosphere 

were obtained from the radio-sonde data, observed at Shionomisaki Meteorological * 


Observatory, about 210km SSW apart from Nagoya, and at Wajima Meteorological 


Observatory, about 240km N apart from Nagoya, at about 0-h and 12-h every day. : 
On the daily means of cosmic-ray in- Gruralation Coefficient Vip tite 
0 =O.2 ~aF ee igs -lLo 


tensities I, barometric pressures P, the 


difference of heights of isobar levels 
< between 100mb and 1000mb Hy, and ps5 
the difference of heights of isobar levels > | b 
: between 100mb and 200mb H,, we take — RS ieee ms ies ee 


me a ' ‘Separe sof Ma shuns 6: ae 5 q 4 . 
res of daily means from their 27-day gf IT Ky | 
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such materials that barometric presstires deviate more than standard deviation of them. 
The coefficients of barometer, effect Brr-.a,7, thus obtained are plotted against the 


corresponding correlation coefficients rrp-#,#, in Fig. 2. From this, such tendency is to 


be seen that the larger the correlation coefficients, 
the larger the coefficients of barometer effect 
become. So we adopt the value 6f coefficient of 
barometer effect corresponding to maximum 


correlation coefficient : 


By=—0.13%/mb (7) 
Bn= — 0.22% /mb Ge) 
By = —0.2223 /mb (W’) 


But, throughout the present analysis, we 
have applied 8y=—0.159%/mb for H, obtained in 
the preliminary analysis at the first stage and 
Brp=—0.24%/mb and By=—0.19%/mb for T and 
W respectively, reported in our previous paper.” 
Hewever, these difference is not so significant 
that the following results on the cosmic-ray 


a0 2 4 6 0 0 mw Ht he Re 


diurnal variations may be changed qualitatively. Local tinas (hewn) 


For example, the semiamplitude of 1st harmonic Fig 3. Diurnal variations of bihourly 
cosmic-ray intensities corrected for 
’ barometer effect averaged over a year. 
averaged in a year, June 1952—May 1953, amounts oe: June 105%ntdy sess 


to 0.7 mb, which gives only the difference of o—o : June 1953—May 1954 
correcting values less than 0.02%, corresponding to the difference of coefficients of 


of solar diurnal variation of barometric pressure, 


barometer effect mentioned above. On the other hand, not so reliable values of 
coefficients of temperature effects a@ra;-paH, and rraz-pa, were obtained. 

3-2. Mean solar diurnal variation of cosmic rays. The solar diurnal variations 
of cosmic-ray intensities averaged over a year are reproduced in Fig. 3, after being 
corrected for barometer effect. In the year, June 1952—May 1953, the total numbers 
of days were 315, 314, and 287 for H, T and W, respectively. In the year, June 1953 
—May 1954, they were 332, 337 and 311 for H, T and W respectively. 

The 1st harmonics of these mean solar diurnal variations are shown on a 
harmonic dial in Fig. 4a. Fig. 4b shows the result of their means taken over those 
two years, giving also that of low energy component (Z) with the range less than 
90cm lead which was obtained by subtracting H from 7. From this result, the phase 
differences between 7, H and L are as follows: eo 

AV(T, H)=30.8+16.6 degrees between 7 and H 
AV(L, T)=41.4+21.1 ms i | eats 
AV (L, H)=72.2 422.0 2 pS ey eehell 

If we were, at the present, allowed to disregard the difference of temperature 
effect between the different components, these phase differences would be attributed to 
the differences of the deflections of primary cosmic rays in the earth’s magnetic field, 
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Fig. 4a 1st harmonics of cosmic-ray Fig. 4b 1st harmonics of cosmic-ray 
diurnal variations are plotted on a diurnal variations averaged over June 
harmonic dial. 1952—May 1954 are plotted on a 


harmonic dial. 


June 1952—May 1953 | June 1953—May 1954 


| | : 
| ce) | O : 
eee 


depending on the mean energies of different components. According to E. Brunberg,» _ 


the mean energy of total components of cosmic rays as giving rise to solar diurnal 


variations amounts to 21-45 Gev, which was measured with narrow angle telescope — 


inclined 30 degrees to the zenithal direction. If we assumes for mean energy of total 
components of cosmic rays T as causing solar diurnal variation to be 21 Gev also, 


though they are measured with narrow angle telescope directed vertically, the above <a 


described phase difference gives the mean evergies of L and H to be 14. Bxt4"4 Gev 3 
and 39.13 Gev, respectively. This mean energy of L will support the discussion in a 


ee 4-2. “ iTS «eG wy Gal eee [sal 
4. Anomalous: solar ira yvaon a eta a somagnetic storms 


a a a a 
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their 24-hour running averages’ and ditirnal harmonics were treated conveniently. 


. Table I 

Numbers of 

Apparatus 
Notation | Location of observatdry . storms adopted 
absorber aperture |70/S1/|S2)| r 
fe Nagoya University = * ee Bias no. +12° Sa LS Re 2om 22 
H * f | 90cm Pb +12° | 34| 15 | 25 | 22 
Ww * a no. +40° | 34 | 16 | 25 | 22 

the Scientific Research “ “ 3 i: 

N3 Institute at Tokyo 72 nD: +85 31/13 21) 21 
NI ; ettee here nel0 cit: Pho} att 0075 ha? 9) 17407) 25 


Taking the day of occurrence of each geomagnetic storms as the Oth day, we 
denote S, for the mean diurnal variation averaged over 3 days before the occurrence 
of geomagnetic storm (from —3rd to —Ist) and S, for the mean diurnal variation 
averaged over 3 days after that (from +1st to +3rd). We difine not S,, but Se Sos, 
as the anomalous diurnal variation associated with geomagnetic storm. Such definition 
was also adopted by S. Yoshida in her analysis. 

4-2. Results. The mean S, averaged over all storms are shown on harmonic 
dial in Fig. 5a.Fig. 5d, to which suffices “ To” are attached. From this it is to be 
seen that vectors S, of different instruments point out 
different directions. In order to confirm whether these 
differences of directions are due to statistical fluctuations 
or not and whether they depend upon the selections: Ese 
group of storms or not, following analysis were caries out. 
At first, N3, having been observed by the group. of 
the Scientific Research Institute is the cosmic-ray intensities 
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with the highest accuracy in our country and 
the statistical fluctuation of its hourly value 
amounts to only 0.26%. So that, we calculated 
Sy of each storm observed by N3, which are 
shown in Fig. 6a. In this figure, dividing 
the plane into two half planes by the straight .. 
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line ZL through original point and perpendicular ‘Fig. 5d Sp of NI Fig 5e. Sp of W 


Fig. 6a Sp of N3 calculated on 

individual geomagnetic storms in 

June 1952—May 1954 are plotted on 

a harmonic dial. The total numbers 

are 37, including 6 of NIJ.. OM is 
mean Sy, vector. 


to the mean S, vector OM, such cases that the end 

points of S, exist in the half plane including the 

point M are 27 (73%) and the other cases are 10 

(27%). That is, we should give attention to the 

following facts as one charactor of S,, “The dis- 

persion of the directions of S, of individual storms 
is very large, but the most of their directions are 
inside --90 degrees from the direction of mean Sp.” 

Next, we computed the mean S, about the 
following groups of geomagnetic storms, on the 
basis of each S, of every storms mentioned above: 

“$1”: dividing into two half planes by the straight 
line L’ through M perpendicular to OM, 
group of such storms that the end points 
of Sp of N3 exist in the half plane not 
including Oz. 

“$2”: among two half planes divided by straight 
line ZL, group of such storms that the end 
points of S, of N3 exist in the half ‘plane 
including M. 


Seite Fi Pena r, for the amplitudes 
ers ae rae Ss ier gh group of 


—aa. -. 
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clear, because diurnal variations ‘observed 
with narrow angle telescopes are large, 
though their accuracies are low: It is 
most clear for N3 because its accyracy is 
the highest.) 

These Sp vectors are reproduced in 
Fig. 7, gathered together into the original 
point. Also, the times of maximum in- 
tensities and the amplitudes of mean S, 


” 


for the group “To” and the group “+r 
are shwon in Fig. 8a and Fig. 8b. From 
Fig. 7, Fig. 8a and Fig. 8b, the following 
facts can be recognized, although the 
errors are large. 


(a) Ampitude. The ratio of ampli- 
tude of S, observed with narrow angle 
telescope directed vertically (ZT) to that 
observed with wide angle telescope (W, 


N3, NI) is about “factor 5.” One reason 


of this is that wide angle telescope measures cosmic 


inlonsitios 4 Splhrun) 


maim 
x 
By 


The time of 
% 


w 


ee ee me 

(¥omPb) = Omay. (degrees) 
Fig. 8a The times of maximum 
intensities of Sp, are plotted 
against the maximum semiangles 


1 
of telescopes. @ BCA WO) ous cater tal 


I 
averaged over all storms. O: 


““y?? = mean averaged over such 
storms that 7a—7»>0 about the 
data of N3. 


Fig. 7 Sp of several kinds of cosmic-ray 
meters calculated on geomagnetic storms in 
June 1952—May 1954 are plotted on a harmonic 
dial. 
O222, ORT wala NS eARANE 


rays incident from different 


Ampltude of Sp (%) 


Sk, sh GO &F “Fo 

(70 cnP6) Oras. (degrees) 
Fig. 8b The amplitudes of Sp are 
plotted against the maximum semi- 


| 
angles of telescopes. e che SET Ores 


mean averaged over all storms. 


d: cigs 


I : 
such storms that 74—7)>0 about » 


the data of N3. 


mean averaged over 
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directions at the same time, so that the amplitude of diurnal variation measured with 
it becomes small for the “average effect.” The other reason is that cosmic rays 
measured with narrow angle telescope directed vertically pass through smaller mass of 
atmosphere and thus contain more contributions from low energy primaries. 

(b) The time of maximum intensities. 

(i) Narrow angle telescope: The time of T is earlier than that of H. 

(ii) Wide angle telescope: The time*of N3 is earlier than that of NJ. 

(iii) Narrow angle telescope and wide angle telescope: The time of N3 is 
earlier than that of T. 

These facts can be explained qualitatively as follws. Let us assume that 
anisotropic flux of cosmic rays exists far outside the earth’s magnetic field and causes 
cosmic-ray diurnal variation as observed on the rotating earth, then the time of 
maximum intensities of diurnal variation depends upon the mean energy of cosmic 
rays observed, owing to their deflections in the earth’s magnetic field. The lower the 
mean energy of cosmic rays, the earlier the time of maximum intensities becomes for 
the primaries of positive charge. Thus (i) is evident. (ii) For N3 is more directional 
than NJ, N3 measures the cosmic rays of lower energies than NJ. (ili) According to 
E. Brunberg’s figures») on the deflection of cosmic rays in the magnetic dipole field, 
WY, which gives the change of longitude of particle’s position is not symmetrical to 
zenithal direction. The difference of ¥ » between direction inclined north and zenithal 
direction is larger than that between zenithal direction and direction inclined south. 
And in the region of low energy, the difference of ¥, between direction inclined east 
and zenithal direction is larger than that between zenithal direction and direction 
inclined west. (For the case of latitude X,,=25 degrees and momentum p=30 Gev/c, 
Y »(48°N)—¥ (0°) =33°, ¥2(0°)—¥ 2(48°S) =5° ; ¥ 2(48°E)—¥ (0°) =60°, ¥ 2(0°)—¥ p(48W°) 
=40.°) Thus, the mean ¥; of cosmic rays measured with wide angle telescope directed 
vertically is larger than that of cosmic rays measured with narrow angle telescope 
directed vertically, so the time of maximum intensities of cosmic- “ray diurnal variation 
measured with the former is earlier than that of the latter. 7 


po. sek a way, the ae: of ic be pees pla Bea of ee teacaret 
< . = J 


ee =: aoa 
the exact « 
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treatment of them is more simple than the others. Fig. 9 shows the differential 
effective primary spectrum calculated by K. Nagashima,” that is, the intensities of 
total components of vertical cosmic rays i(E, x,) produced by the primary cosmic rays 
of kinetic energy E at the atmospheric depth x,=10mH,O. Curve I and II correspond 
to two limiting cases. In the same figure, the deflections of primary cosmic rays in 
the earth’s dipole field are reproduced from E. Brunberg’s figures.» For T, we take 
magnetic cut-off energy in Japan, geomagnetic latitude 25°N, as 10Gev. Next, we 
must determine the effective primary spectrum of high energy cosmic rays H. The 
effective primary spectrum shown in Fig. 9 was computed on the ground of latitude 
effect only and there is not significant difference of latitude effect of vertical cosmic 
rays between total component and hard component.” Therefore, we assume that the 
shape of effective primary spectrum 
of hard component is the same as 
that of total component shown in 
Fig. 9, with the exception of the 
difference of absolute flux. Taking 
the intensity of T as 100%, the 
intensities of hard component and 
H are 70% and 50%, respectively. 
Thus, 20/70=29% of hard com- 
ponent are absorbed by the 90cm 
lead. So that, in Fig. 9, we draw 
the inclined straight line from the 


; : Fig. 9 Differential effective primary spectrum 7(E, x.) 
point of Z=10Gev in such a way at sea levél is plotted against kinetic energy E of 


that the low energy part of curve primary proton. Curve I and II correspond to two 
limiting cases. (after K. Nagashima) The change of 
longitude of particle’s position by the deflection in the 
of total intensity. Now we can earth’s magnetic field is reproduced in the same figure. 


consider the high energy part as (after &. Brunberg) 
the effective primary spectrum of H. 

Here, according to T.H. Johnson,” it is known that latitude effect and east-west 
effect of cosmic rays vanish at the depth of 16mH,O, namely, field sensitive rays of 
primary energies from 7 Gev to 14 Gev are absorbed at that depth. On the other 
hand, 90cm lead is just equivalent to 6mH,O of air for the range of ~ meson. Thus, 
almost all of secondary particles produced by primary cosmic rays of energies from 
10 Gev (cut-off energy) to 14 Gev are absorbed by 90cm lead. This is the reason why 
we draw the inclined straight line in such a way that the part of low energy is 
absorbed in order to determine the effective primary spectrum of H. Also, we have 
seen in subsection 3-2 that the mean energy of low energy component L defined as 
the difference between Z and H is about 14.5 Gev. This result supports above 


consideration too. Paes 
Let us assume that anisotropic flux of cosmic rays is formed far outside the 


earth’s magnetic field by any mechanism of acceleration in a certain direction, by — 


cut off by this line amounts to 29% 


ety a) ual aoe Reaiiari 
xh 
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which cosmic rays obtain a gain of energy 4EGev. Then, as K. Nagashima® pointed 
out at the first time, the change of intensities are: 


2 : 1+4ak* 2(E+m,) 


] a 0 1(E, ’ 1 
AE x) 4E-{ ee ee er i(E, x2) (1) 


a=00.9, constant number 
=0,93 Gev, rest energy of proton, 


where the terms of orders higher than (4£)? are neglected and the primary spectrum 
obtained by H.V. Neher? was used. The first term in bracket arises from the shift of 
energy spectrum and the second term arises in connection with Liouville’s theorem. 

Taking the space distribution of anisotropic flux as A,=4i(E, 0) cos(¥) and 
denoting the change of longitude of particle’s position by the deflection in the earth’s 
magnetic field as ¥,(E), diurnal variation on the earth is: 


A=| ” ail, x,)cos(¥ —¥ (E)}dE/\” i(E, x,)dE 
=A(E,)cos{V —¥(E,)} ~ 
AXE.) =| \" 4i(E, %,)cos v (EE } +f (" 4ite, x,)sin¥ (E)dE } /\ ie «ee 
fan’ EJS ie Ai(E, x,)sin¥ (E)dE ie AE, X,) cost E)AE (2) 


Taking he 6x10-* Gev, the results caculated from eq. (1) and eae (2) are given rf 
Table II. It gives 24~26 degrees as the difference of phase of diurnal variation ; 
between T and H. On the other hand, from Fig. 7 and Fig. 8a, the phase differen 
of Sp between T and H is about 60~90 degrees, which is considerably larger than the 
coset ihe Rg a the experimental error is ws 
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on a harmonic dial in Fig. 10a. ” Frequency distribution of the time of maximum 
intensities of those is shown in Fig. 10b. The time of maximum intensities of mean 
Sp is 9.5 hour in local time. 


The Leal Leet of retina intimate of Cobmmic any dinrmal variation T(heuss) 
te cs i. a ae 


= Bs sf 


et fs 


Chet ei pa es ee ee 


Oo 310 4 40 fe ie me 760 fe see ate ig ile al¢ a0 sae Whe Sho @ 


The phase of macimatm intnaitica of rtmic-ray diurnal variation ¢ (ugha) 
Fig. 10b Frequency distribution of 
the time of maximum intensities of 


Fig. 10a Each Sp of N3 calculated on such each Sp of N3 calculated on such 
geomagnetic storms in July 1947—May 1954 geomagnetic storms that 7a—7>0. 
that 72—75>0 are plotted on a harmonic dial. =—-—— -> Mean 


OM is mean Sp vector. 


——-—_ ~~" = 


Next, we calculated S, in every year averaged on such geomagnetic storms that 
the amplitudes of diurnal variations increase and the ranges 4H of horizontal force 
exceed 1007 respectively. They are shown in Fig. lla and Fig. 11b. The numbers 
of geomagnetic storms in each case are given in Table III. From these, it is to 


—— 5 


ee 


N 
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seen that S, in every year point out about a definite direction and shift their position 
(S,) from year to year. The movement corresponds to the change of phase of cosmic- 
ray diurnal variation with 22-year’s cycle found by T. Thambiyahpillai and H. Elliot.1% 
Thus S;, namely, S, after S. Yoshida” is considered to consist mainly of component of 


secular variation on the average. 


Table, U1 

ee £. i. Clacsiivavon The une of geomagnetic storms 

Period it = peat AH>100y 1s —Ta>0 Total : 

Jul., 1947—Jun., 1948 1 19 ok j 

Jul., 1948—Jun., 1948 26 36 
Jan., 1950—Dec., 1950 24 18 
Jan., 1951—Dec., 1951 23 24 
Jan., 1952—Dec., 1952 21 20 
Jan., 1953—Dec.,, 1953 - 13 12 
Jan., 1954—May., 1954 0 10 
77 months 119 130 


~ This result leads us to the following physical picture. Anisotropy S, is already 
formed in a large space around the sun in a long range from the earth. At the time > 
of geomagnetic storm, the earth is immersed in a perturving medium, thrown out from 
the sun and generating geomagnetic storm. Then cosmic rays are accelerated when 
passing through the perturving medium, and anisotropy S, is formed in the neighbour- 
hood of the earth. So we observe the resultant of S, and S, at the time of 
geomagnetic. storm. a 


i; 6. Long-time change of diurnal variation of cosmic rays 


12 months’ averages of first harmonics of cosmic- -ray diurnal variations were — 
computed in every month, using the data of T and W from March 1949 to Febhiig a 
1956 ae the data of V3 from January 1950 to August 1955. (Among these series, the 
observation from October 1951 to May 1952, that of W from October 1951 
: July 49 were 
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Fig. 12a The amplitudes of 12 months’ averages of Ist harmonics of cosmic-ray 
diurnal variations are plotted in every month. 
e—e: T, o—o: W, x—x: N3 
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Fig. 13 12 months’ averages of K-index K (monthly mean of daily sum) and that of 
relative sunspot number R (monthly mean) are plotted in every month. 


7. Summary 


1) Mean solar diurnal variations of cosmic-ray intensities of total component — 
T and high energy component H penetrating 90cm lead absorber measured with narrow 
angle telescope directed vertically show the phase difference of about 30 degrees from 
each others. 
2) Though it is well known that Seomaloas solar diurnal variation Sp of cosmic 
rays appears at the time of geomagnetic storms, the existence of Sp about T and HA 
were confirmed also. i. 
st , 3) Amplitude: of ce of T is the largest and about five times as eS as 
measured with wide angle icles iy : tes i, ) 
5 ie aiks, Oca ee 


i 
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1952 and its increase is especially remarkable for 7. It may be connected with the 
remarkable increase of K-index, not being rather asociated with relative sunspot 
number. 

We intend to continue this experiment with narrow angle telescope directed 
vertically with higher accuracy pial better stability than the present one. 
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Abstract 


An experimental method is described, designed chiefly to furnish data for 
a study of the direction, and stability with time, of thermoremanence in rocks. 
Specimens were heated to the Curie point in an evacuated electric furnace. Two 
pick-up coils were arranged close to the gap of a tuned a.c. electromagnet providing 
a maximum field H of 2400 oersteds. These were balanced to make their 
resultant e.m.f. zero in the presence of H alone, and proportional to d//dt when 
a specimen of intensity of magnetization J was in the gap. This e.m.f. was 
applied to the vertical plates of a cathode ray oscilloscope. The potential drop 
over a small resistance in the electromagnet input was applied to the deflection 
coils, giving a measure of H. Computations based on the resulting pattern on 
the c.r.o. screen yield a loop of the /-H type, with J in arbitrary units, from 
: which the coercivity can be evaluated. A pronounced ‘‘ sawtooth’’ pattern has 
been observed in the (dI/dt)-H traces of pyrrhotite and franklinite specimens, 

particularly just below the respective Curie points. 8 


a mae ~The Continuous Recording of Hysteresis ¢ Curves 


Pe investigations [1], [2], [3], [4] have been made into the magnetic bytes 
of rocks, gees by stepwise measurements using a ballistic arran t 


eTeiaand 
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integration of the vertical sweep was required to convert this into an --H loop. 
Igneous rocks were investigated.at room temperature in maximum fields of 750 


oersteds. 
Experimental Requirements 


It can be shown that the fiéld required to saturate spherical grains of magnetite 
falls in the range of 2000 to 3000 oersteds at room temperature. For several other 
strongly ferromagnetic minerals the saturation field appears to be also of this order. 
The equipment should therefore be designed to provide a magnetizing field of 2000 to 
3000 oersteds. This field, as well as the temperature, must be as nearly possible 
constant over the volume occupied by a standard specimen, ie. 4 to 5 cm’ in the 
equipment described here. 

The method must be capable of measuring the relatively small intensity of 
magnetization of a specimen, in the presence of the magnetizing field. For a standard 
specimen of incremental volume susceptibility 2000x10-*c.g.s. units, measured in a 
weak field, this intensity accounted for roughly one part in 2000 of the total induction 
in the equipment employed. 

To make observations up to the Curie points of the materials concerned, a 
maximum temperature of 600°C was necessary for the ferromagnetic minerals here 
examined, and this involved a satisfactory method of measuring the temperature. To 
prevent oxidation of the specimen in the laboratory, it must either be surrounded by 
an inert atmosphere, such as nitrogen or argon, or heated in a vacuum. 


Summary of the Method 


The method used by Bruckshaw and Rao, briefly described above, was modified 
for work at high temperatures. A tuned circuit at 50c/s in which an electromagnet 
acted as the inductance now provided magnetizing fields up to 2400 oersteds. Two 
coils were placed near the air gap, the first being mounted on the core. The second 
coil, in an adjustable position outside the core, was threaded by part of the flux 
spreading outwards from the gap. These coils were balanced in such a way that the 
e.m.f.s set up in them by the magnetizing field alone, cancelled out. 

The specimen was in a vacuum inside an electric furnace. With specimen and 
furnace in the gap, the alternating field magnetized the specimen, causing a negligible 
increase in the total flux of the magnetic circuit. Due to the respective positions of 
the pick-up coils in relation to the gap, a differential e.m.f. was produced, and was 
directly proportional to d//dt of the specimen. 

This e.m.f., after suitable amplification, was applied to the vertical deflection 
plates of an oscilloscope. The potential drop along a small resistance in series with 
the electromagnet was applied to the horizontal deflection coils of the oscilloscope, 
giving a measure of H. os 

Again it was necessary to convert the vertical scale from dJ/dt to J, and 
these calculations were made on data obtained from an enlarged photograph of the 
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pattern, yielding a hysteresis loop 
of the LH type, in which J was in 


arbitrary units. 
The Magnetizing Field 


Electromagnet (Fig. 1) 

To avoid eddy currents, the core 
was built up from suitably trimmed 
transformer laminations of “ stalloy,” 
a high-permeability silicon steel. A 
3cm gap was left between the pole 
pieces, which were interleaved with the 
yoke to prevent breaks in the magnetic 
circuit. 

Two magnetizing coils VW, each of 
1000 turns of no. 16 S.W.G. specially in- 


DruTscH 


cm Z 
CMs io “|(b) 


Fig. 1 


sulated copper wire, were wound on Tufnol formers of rectangular cross-section to fit the 


yoke in the position indicated in Fig. 1. 


The complete electromagnet was assembled 


with Tufnol clamps and bolted to a wooden table in a horizontal position. 


Tuning Circuit (Fig. 2) 


Amp 


The power for the magnetizing 
coils was obtained from the mains 
through an auto-transformer which 
supplied a maximum of 270 V, insuffi- 
cient to drive more than a small frac- 
tion of the required current through 
the coils of total 
The 
condenser C in series with the coils. 
Since the permeability of the core varied 
with the field, the self-inductance was 
a function of the magnetizing current, 


impedance 720.2. 
impedance was reduced by a 


and a range of capacities was necessary 
to obtain different peak fields. 

Since the gap was wide, the 
largest flux-density attained by the 
electromagnet was only 1700 G, at 
3amp (both r.m.s. values), corresponding 
to.a field of amplitude 2400 oersteds. 
All fields up to this maximum could 
be obtained by tuning the magnetizing 


coils with the range of capacities 3.7 


a 
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to 43 nF. 
Measurement of the Magnetizing: Field 

A search coil whose effective area and shape were approximately those of a 
standard specimen could be positioned centrally in the gap. The reading on an a.c. 
voltmeter, multiplied by the appropriate factor, then gave the average field amplitude 
Hy over the specimen volume. ~ - 

A gap survey, carried out with a smaller search coil, revealed that the field in 
the volume occupied by a standard specimen varied by 42% from this average, the 


maximum field being found along the line joining the centres of the pole pieces. 


The Pick-up Circuit 


A rectangular former F wound with 200 turns of no. 30 S.W.G. enamelled copper 
wire was fixed rigidly to one of the pole’pieces, the coil centre being 3.7cm from the 
gap centre. 

The second, adjustable pick-up coil A, with 2500 circular turns of double silk- 
covered wire (40S.W.G.) was placed close to the gap, to be threaded by some of the 
flux spreading outwards. As the ratio of the flux lines threading each of the two 
coils varied slightly for different magnetizing fields, the air-cored coil could be 
balanced by rotation with a fine adjustment control. Both formers and the adjusting 
mechanism were made of Tufnol. 

Harmonics Coils 

When the fundamental frequency was balanced a small e.m.f. remained, consist- 
ing at the lower fields chiefly of third harmonic. The relative positions of the two 
pick-up coils (Fig. 1) were chosen by trial to give a minimum of these harmonics. 

A large coil H1 of 2100 turns was mounted in a position (Fig. 1) where the 
resultant flux due to the field surrounding the two magnetizing coils was zero at the 
fundamental frequency. The e.m.f. set up in this coil then consisted mainly of third 
harmonic. After suitable adjustment of the amplitude by a potentiometer in the 
circuit, this was balanced against the harmonic residue in the pick-up coils, most of 
which was eliminated by this method. In fields above 1500 oersteds, a small coil H2 
of 900 turns (Fig. 2) reduced the amplitude of the undesirable harmonics, but since 
these were now more complex, they balanced only to a minor degree. Both harmonics 
coils were made of double silk-covered copper wire (40 S.W.G.) wound on Tufnol 
formers whose position above the electromagnet was adjustable. 

All the leads of the pick-up circuit (Fig. 2) were electrostatically screened. 


The Oscilloscope Pattern 


The oscilloscope was a Cossor double beam instrument giving a blue screen 


pattern. 


Vertical Deflection. atebze 
The output of the pick-up system was fed into the c.r.o. amplifier, with a 


maximum gain of 900, and the amplified e,m.f. applied to the vertical plates of the 
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tube. 


Horizontal Deflection. 
The potential drop across a small resistance FR in series with the electromagnet 


was applied to the c.r.o. deflection coils D. The plates could not be used since both 
pairs, as well as the autotransformer, were earthed on one side. Depending on the 
input current, the tapping was made from one of a set of nichrome resistances in the 
range 0.3 to 1.72, selected to give full screen deflection. 

By means of suitable switches x and y in the c.r.o. input, the coordinate axes 
could be superimposed on the photograph of each pattern. 

Deflection Calibration. 

The deflections on the tube screen were found to be proportional to the input, 
up to the dimensions of the pattern. These were generally confined within a central 
area of 5cm diameter, where the screen was made flat to eliminate curvature error. 
Photography 

A Cossor 35mm camera with manually operated shutter was used. Blue-sensitive 
Kodak R20 and Ilford 5BIl films were selected, both types having high resolution. 
With low brilliance to render fine outlines on the curves, a suitable exposure for 


2 s. The axes required much shorter timing. 


patterns was 7 


Furnace and Thermocouples Chrome! 


A 
The furnace (Fig. 3), essentially a Lome 


long quartz tube wound with resistance 


wire, was of the general type used by 
itri—— Furnace top 


Tomlinson and Bockris [11] in the range (pyrex) 


1000 to 2000°C. It was sufficiently narrow 
to fit between the pole pieces, leaving an 


Thermocouple 


air cushion 2 to 3mm wide on each side. element 
A short, perforated quartz tube, Voce = ie beets 
quartz 


fused inside to the bottom of the furnace, outlet 
supported a quartz crucible containing Mgin quartz 


Positioning rings 
the specimen. Heat conduction away tube 


(asbestos) 
from the crucible was thus greatly reduced, To a.c.. 
whilst convection was overcome by supply 
evacuating the tube through an outlet. A 


ground pyrex cone sealed the furnace Crucible 
containing 


vacuum-tight at one end, and, since the ‘ 
Specimen 


Crucible 
support 


vacuum grease melted after slight heating, 
the furnace top was water cooled. 
Tungsten rods sealed through the pyrex cm 
provided contacts for the thermocouple. b—"— 
The furnace was wound with 43 ° s 10 
turns of no, 27 S.W.G, “nonmagnetic ” Fig. 3 


wae 


The Measurement of Magnetic Hysteresis in Rocks and Minerals at High Temperatures 113 


nichrome wire, concentrated at the top and bottom. With this arrangement, the 
temperature in the region normally occupied by a specimen was found to vary by 
+1% from its average value. A sheet of mica surrounding the Wiring section, together 
with the air cushion, provided partial heat insulation for the pole pieces. Since, at 
high temperatures, the tube was left in the gap for at most 15s., further insulation 
was found unnecessary. The a.c. supply from an auto-transformer heated a specimen 
to 600°C in 10 min, with 2.75 amp, but no advantage was gained from heating so 
quickly. To prevent induction in the pick-up coils from the heating coil, the supply 
was cut off for a short interval during the photographic exposure. During this interval 
the specimen temperature dropped by about 1°C at 600°C, and this was neglected. 

The furnace was kept in a raised position until the required temperature was 
reached, and then lowered into the gap by means of a brass ratchet arrangement. 
Thermocouple Circuit 

A chromel-alumel element yielding a large thermal emf. was used. The 
junction was arc-welded into a bead, in a reducing atmosphere. To diminish the 
magnetic effect caused below 140°C by the nickel in chromel and alumel, the 
thermocouples were made of no. 26 S.W.G. wire, the smallest gauge available at the 
time. Any remaining magnetic background was accounted for by subtraction of the 
oscillogram pattern due to the thermocouple alone, from that obtained jointly with 
specimen and thermocouple. This precaution was only necessary with a few relatively 
weakly magnetic specimens. The hot junction was positioned centrally inside the 
furnace (Fig. 3). 

A chromel-alumel cold junction was placed in ice inside a Dewar flask. Readings 
were taken on a galvanometer whose full deflection corresponded to about 60°C, the 
higher temperature ranges being obtained by using suitable resistances in the circuit. 
To prevent secondary thermal e.m.f.s., all the wire in the thermocouple circuit was 
thick chromel or alumel. Pairs of metal contacts at the galvanometer, furnace top, 
cold junction head and switches, respectively, were closely spaced to equalize the 
temperature of each pair. 

Accuracy of Temperature Measurements 

A standard platinum-platinum/rhodium couple was used to calibrate two chromel- 
alumel elements. The latter gave readings differing by up to 2°C at 600°C, where the 
temperature could be read to 1°C. Together with smaller errors, due to local heating 
at contacts, lag in the galvanometer response, etc., the uncertainty in the temperature 


determination may have amounted to +5°C at 600°C. 


The Vacuum System 


The furnace was evacuated by means of a two-stage pump, the vacuum circuit 
being thick-walled rubber tubing sealed to glass connections with wax. The vacuum, 


measured with a mercury level gauge, averaged about 1mm of mercury, sufficient to 


prevent any significant oxidation of specimen or thermocouple in the temperature 


range used, 
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The Specimens 


Specimens investigated included the strongly ferromagnetic minerals magnetite, 
pyrrhotite and franklinite, as well as rocks such as basalts, in which magnetite or 
titanomagnetite commonly occur as minor constituents. Specimens were cut into 
prisms 3.0cm long and 1.1x1.1cm in cross-section. When powder concentrations of 
a ferromagnetic mineral were required for study, the mineral powder was mixed with 
plaster of Paris powder. Water was added and the resulting paste allowed to set in 
a mould. Artificial “rock” specimens were thus produced. 


Computation of the Hysteresis Curves 
If H denotes the field in the magnetic circuit and B the induction, then H=B 
in the air gap. A specimen sets up further lines of induction, such that 
B=H+4a7l @® 
where / is its intensity of magnetization. Since B, H and J are vectors, this relation 


is numerically exact only for specimens saturated in a uniform field. Using an a.c. 
electromagnet, the e.m.f. set up in a pick-up coil is then proportional to dB/dt, where 


dB/dt=dH/dt+4n(dI/dt). ~ (2) 


Since dH/dt is balanced out in the method used, the field term vanishes and the net 
e.m.f. becomes proportional to dJ/dt only. For a sinusoidal input current of frequency 


fi 
H=Fhsinot pret 


where H, is the peak field and w=2zf. 
So 
dH/dt=H,wcosot, (4) 
and since esd 
~ _ al/dt=(dI/dH) x (dH/dt) jonktbebiid 


dlI/dte (dI/dH) cosat. ‘yp Temivenie ee ihe 


An ideal pattern on the c.r.o. screen is then proportional 


ponenely and to (6) vertically. In practice, the. length | of the orzo 
inc H up to nearly 1200 Hel: yaaa Pegi - fh ie 
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oscillogram, taken at H,=2300* oersteds, has been retraced in Fig. 4. At each 
interval, the difference between the main loop (a) and loop corresponding to unbalanced 
harmonics (b) was measured, giving net readings caused by the specimen alone. Values 
proportional to d//dt thus obtained along the ordinate were converted into quantities 
proportional to d//dH after division by cost (Fig. 5). A numerical integration between 
limits of 0.05, proceeding from “the origin at an arbitrary abscissa, then yielded the 
j desired /-H curve (Fig. 6a). 

Half the vertical separation between 
the maximum and minimum of J was taken 


as the saturation magnetization I,. Half the 


separation between the intercepts at H=0 


Fig. 5 


gave the remanence I,. The horizontal axis 
could be obtained by locating a point along the ordinate axis either halfway between 
+I, and —J,, or halfway between +J, and —J,. In an ideal hysteresis loop these 
points should coincide at the true origin, but due to experimental errors involving 
mainly the oscilloscope measure- 
a REE ments, the loop may be slightly 
distorted and fail to close. When 
i this has been the case (Fig. 6b), 


the final abscissa was drawn 


21 : through an origin located in a 
s 


-2000 -I000 
t+ 


mean position along the ordinate. 
e. plier Se: Half the separation between the 
intercepts on’ the abscissa then 
Nl Ay yielded the coercivity H.. 


From the distortion in the 
Fig. 6a 
quoted example, H,, J, and J, 
were estimated to be accurate to +1.6, 2.5 and 0.6% respectively. These may be 
regarded as typical figures for strongly magnetic rocks. In smaller fields, involving 
simpler patterns, these errors decrease, 
whilst measurements at high temper- 
atures or with weaker rocks entail 
correspondingly larger errors. Just 
below the Curie point, the /H loop 
tends to become very flat, and then 
the expected error in measuring H, 
may be as large as 10 to 20%. 
The “Sawtooth” Pattern 
In Fig. 7 (a—f), oscillograms are 
reproduced for a pyrrhotite specimen 
(P1) at H)=1200 oersteds and temper- 
atures between 20°C and 310°C. To 
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facilitate the computing, the harmonics loop was expanded to surround the main 
pattern, which remained relatively unchanged up to 210°C, but rapidly altered between 
295°C and 310°C, corresponding to the collapse of the /-H loop just below the Curie point. 

At high temperatures, the oscillograms exhibited the “sawtooth” pattern (Fig. 7 
(c—f)), which was most pronounced in the region of the coercive field and affected 
those parts of the loop where magnetization changes are likely to proceed chiefly 
through irreversible domain movements. Thé'pattern resembled a damped oscillation : 
with a regular period. 

A phenomenon possibly analogous to this has been observed by Cross [12] in the | 


neighbourhood of the electrical coercive 


field of barium titanate crystals, where 
a “kinking” occurred at a certain 
temperature below the Curie point. 
In the present case, however, the 
“sawtooth” pattern was not noticeable 
in the computed LH loops, except 
when they were greatly enlarged. 

Three other examples of the 
“sawtooth” pattern, retraced in Fig. 
8 (a—c), appeared to be most pro- 
nounced in the range 15°C to 50°C 
below the Curie points concerned 
(Table I). 

Results of hysteresis measure- 
ments on a number of specimens are 
discussed in a separate paper in this 
issue of the Journal. 


ons 


how 


“phy i eee ye ese iy 
. >. ; aie x 

o tf peat be OP Pes) 
= “ON eae 


The Measurement of Magnetic Hysteresis in Rocks and Minerals at High Temperatures 117 


Table 1. Temperature’ of pronounced “ Sawtooth”? patterns 


Specimen ah wbetcC) ES) Die) 
P1 (Pyrrhotite) 330 305- 300 25-30 
a 
P2 (Pyrrhotite) ; 320 295-270 25-50 
F 4 (Franklinite) j 20 5 15 
— 
F5 (Franklinite) 525 480 45 


T.=Curie point. 
T;=Temperature(s) at which pattern was photographed. 
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Abstract 


The coercivity, remanence and saturation magnetization of a number of 
minerals and rocks have been obtained as a function of temperature, using 
alternating fields up to 2300 oersteds. These results suggest that most of the 
specimens examined can retain thermoremanence for long periods, even at high 
temperatures. Irregularities occurring in some of the temperature curves indi- 
cated the presence of more than one ferromagnetic constituent. Heat treatment 
in several cases caused considerable changes in the hysteresis curves measured 
at atmospheric temperatures. The thermomagnetic behaviour of two basalts has 
supported the view that the reversed magnetization of certain rock formations 
is due, not to abnormal magnetic properties on the part of the material concerned, 
but to a past reversal of the geomagnetic field. 


Introduction 


As a result of field and laboratory work undertaken since 1946 by the Geophysics 
Department of the Imperial College, it has been discovered that the Tertiary tholeiite 
dyke system of Northern England is permanently magnetized in a direction opposite 


that of the present earth’s field [1]. It was subsequently shown [2], [3] that many of the. 
Tertiary lava flows and dykes of the Isle of Mull possess statistically significant mean _ 


directions of magnetization, some of which are in opposition to the geomagnetic field, 
whilst others nearly coincide with it. Reversals in the polarization of igneous rock | 
formations have been discovered in many countries, including South Africa, Iceland, 
Japan, France, etc. ae 

Igneous rocks usually acquired their oS residual magnetism while coc 
in tS earth's gee field from | temperat OV 


a 
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As an alternative to this hypothesis, Néel [4] and others have suggested that an 
intrinsic property of the rocks themselves may account for their reversed polarization. 
In support of this, a number of Japanese specimens do acquire thermoremanence in 
the laboratory, in opposition to a small applied fielc 1[5], [6]. An ultimate decision in favour 
of one or the other of the two basic theories thus requires further research into the 
magnetic and thermomagnetic properties of rocks 

In a paper by Bruckshaw and Rao, [7] experiments are described in which the 
magnetic hysteresis of rock specimens, mainly from the tholeiite dykes, was measured 
in fields of 750Oe. Large coercivities were obtained, comparable to those of certain 
magnet steels. These data permitted the assumption, essential to the field reversal 
hypothesis, that the inversely magnetized Tertiary rocks have retained a large pro- 
portion of their original polarization for 30 million years, perhaps mainly in opposition 
to the earth’s field. 

The present paper deals with the results of an extension of the above work 
into a study of the magnetic hysteresis of minerals and rocks at temperatures up to 
their Curie points. The experimental method used is the subject of a separate paper 
in this issue of the Journal. The knowledge thus gained can allow deductions to be made 
about the stability of thermoremanence. For example, the coercivity of a basalt 
with its Curie point at 580°C may exceed —100 oersteds in the entire range from 
500°C to atmospheric temperature. It is then safe to conclude that the rock will be 
capable of retaining the direction of its original thermoremanence not only for long 
periods and in opposition to the ambient field, but even if partial reheating should 
subsequently occur. Very few data on this subject have been published. Koenigs- 
berger [8] measured the coercivity of a hematite specimen up to 500°C and found that 
it decreased rapidly at first, but more gradually at the higher temperatures. Akimoto [9] 
reported a similar behaviour for the case of the ferromagnetic constituents separated 
from an andesite sample, while Forrer [10] observed that the coercivity of a magnetite 


crystal decreased steadily up to its Curie point. 


Results 


If the specimens were saturated at room temperature in a field of 2300Oe., the 
usual terms coercivity (H,), remanence (I,) and saturation magnetization (I;) will be 
used. If saturation was not achieved at 20°C, the corresponding terms coercive force 
(H,!), residual magnetization (I,') and maximum magnetization (I,’) will be used at 
all temperatures. iz 

The rock and mineral samples were cut into cectapaniay prisms of approximate 
1.2cm2 section by 3.0cm length, the field being applied perpendicular to the long axis. 
Magnetite, pyrrhotite, franklinite and one rock type, basalt, were investigated. Except 
for one powdered magnetite, all specimens were cut from natural rock samples. 

(1) ~Magnetite 

(a) Magnetite from Pretoria, South Africa, (specimen M1) was ORenIns at 20°C 

and then at a sequence of temperatures as the rock cooled from 604°C. The maximum 
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field H, was 23000e. in all cases. 
The results (Fig. 1) show the Curie point of the specimen in the range 580-585°C, 


in good agreement with the value for pure magnetite. The saturation magnetization 
shows its sharpest drop between 400 and 435°C. On the basis of the domain theory, 
I, represents the sum of the spontaneous magnetizations of all the domains in a 
material and, for pure magnetite, the curve.of-J, against 7 should be smooth, J, 
falling more and more rapidly as the Curie point is approached. Specimen M1 there- 
fore is likely to contain a second constituent with its Curie point in the neighbourhood 
of 435°C, possibly a titanomagnetite, which is the name generally given to a magnetite 
containing TiFe,O, (ulvéspinel) in solid solution. Since ulvéspinel is paramagnetic, its 
concentration in the solid solution determines the magnetic properties, including the 
Curie point, of the titanomagnetite. 

The temperature curves were all irreversible (Fig. 1), the value of J, at the end 


of the experiment being only 17% of that determined before the specimen had been 
heated. Akimoto [11],[12] also obtained irreversible 7-T curves. These results suggest 
that specimen M1 initially contained maghemite (y-ferric oxide) which turned into 
the o-form (hematite) after heating. Due to defects in its crystal lattice, maghemite 
can be strongly ferrimagnetic, but it becomes unstable when heated and turns into 
hematite at 275°C, according to some observers, and at 400-800°C, according to others. 
The Curie point of both ferric oxides is given as 675°C. Hematite is possibly 
ferromagnetic, with low values of J, and J, but high values of H., possibly several — 
thousand oersteds [13]. This would account for the uncommonly large coercivity (—613 Oe.) 
exhibited by the specimen at 20°C after heating. Koenigsberger [14] also examined 
magnetite from Pretoria and found H,= —380 Oe. for a solid sample and H,=—450 to 
—570 Oe. for a powder. On the other hand, he observed that “titanomaghemites % 
yielded principally magnetite rather 
than hematite after prolonged 
heating to 700°C. The fact that no 
Curie point above 585°C was ob- 
| served in specimen M1 suggests 
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and after heating to 600°C on fotir successive occasions. By far the greatest changes 
resulted from the first heat treatment (Table I), and these were of the same order of 
Poanie Yo Sneed AT 2 en aneee We tok magnitude as in specimen 1/1. 


treatment : (c) As several observers have 

T=20°C. N=Number of times the specimen 
has been heated to 600°C ; 
—neneene ee reer ar A eh POcKsS Cancbe ‘reduced essentially to tithes 


pointed out, [15] the magnetic properties of 


oss] oersteds arbitrary units of the ferromagnetic minerals contained 
| N ‘AS | ik: Te in them. Magnetite from Tintagel (Corn- 
Yt se Ae | e a eS 5 <8 
| 0 | 973 ] 46.0 153.1 | wall) was powdered and the magnetic 
Votie| 52 1623 Chia! constituent extracted with a magnet. 
bt? bet G72 8.2 10.6 Sufficient of the ungraded powder was 
3 | 682 7.9 | 9.3 added to plaster of Paris to produce an 
| 4 | -607 8.0 14.0 A sl 

artificial specimen containing about 25% 


by weight of magnetite. Measurements 
were made at temperatures rising from 20 to 600°C, with H)=2300 Oe. 

The Curie point was again in the neighbourhood of 585°C (Fig. 2), but the J,/-T 
curve here resembles the /,-T curves of a pure ferromagnetic material. The H,’/-T 
and J/,/-T curves exhibit a peculiar feature, both showing a rapid rise up to about 
100°C. This effect may well be due to incomplete saturation at normal temperature. 
The location and magnitude of the resulting maximum would then depend on the 
manner in which /,’, J,’ and H,’ vary with field and temperature, respectively. 

(2) Pyrrhotite 

(a) Specimen P1, of undetermined origin, was examined upon cooling from 
345°C in a field H)=12000e. With increasing T (Fig. 3), the decrease of H,’ is more 
gradual than that of either J,’ or J,’,-but-all three curves fall to zero at 333°C. The 
sharpest drop in J,’ and J,’ occurs near 300°C, so that the region between 300 and 
333°C possibly includes two Curie points. The LH loops in Fig. 4 (a—d) show the 
rapid increase of hysteresis for a small temperature drop just below 333°C. 


25 25 
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Fig. 4 Specimen Pl. Pyrrhotite. 1-H, Curves. 


At 20°C the hysteresis loop showed an approach to saturation. “Several studies [16] 
have revealed the strong anisotropy of pyrrhotite crystals, which in certain cases 
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attain saturation only in very high fields. Davis [17] investigated powdered pyrrhotite 
which was not saturated at 33400e. On the other hand, the hysteresis loop of a 
pyrrhotite specimen obtained by Tzu-Chang-Wang [18] shows near-saturation at 501 Oe., 
H, being of the same order as in specimen P1. 

(b) Specimen P2, from Sudbury (Ontario), retained ferromagnetic properties 
beyond the Curie point of pure pyrrhotite and was examined as it cooled from 590°C 
in a field H)=2300 Oe. (Fig. 5). The J,/-T curve terminates at 560-580°C, which is 
near the Curie point of magnetite. Its shape suggests that this curve contains two 
“simple” components of magnetite and pyrrhotite, respectively. Below the Curie 
point of pyrrhotite the curve is thus a resultant of the J,/-T or J,-T curves of each 
of these constituents, whilst above that temperature it represents magnetite alone. 
From the intersection of tangents drawn to the J,/-T curve, the Curie point of the 
pyrrhotite was estimated to be 320-325°C. The irregular trend of H,’ near 300°C may 
‘have arisen as follows :— 

At room temperature, the measured value of H,’ is intermediate between the 
separate coercive forces of pyrrhotite and magnetite. A comparison of Figs. 2 and 3 
suggests that, in the neighbourhood of 300°C, H,’ of the magnetite is larger than that 
of the pyrrhotite, but that they both decrease with rising temperature. However, the 
contribution of the pyrrhotite to the total magnetization also diminishes until, at its 
Curie point, the value of the resultant coercive force is a maximum and represents 


magnetite alone. It follows that a minimum exists below this temperature, as confirmed — 


in Fig. 5. A similar mechanism pote account for the observed irregularity in the 
L,/-T curve. 
(3) Franklinite 
All franklinite specimens have come from the Franklin he Jersey) area, and 
Fl to F4 have been cut from a single sample. 

_ (a) Specimen F'1 was investi- | 
gated as it cooled from 605°C in 7% 
a field H)=12000e. The graphsin _ 
Fig. 6 suggest a Curie point just =, _ 
above 600°C. The H,! maximum “p> 
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between 61 and 20°C, where J,’ increased by a factor of 4.4 and H,/ dropped sharply. 
Since ten hours had elapsed between the measurements at 61 and 20°C, respectively, 
it was considered possible that these marked changes in magnetic properties might 
depend on time. 

(b) Specimen F2 was investigated during the heating as well as the cooling 
process. The field was increasefl to 2300 Oe, sufficient to cause saturation not far 
above room temperature; hence the location of both coercive force maxima at relatively 
low temperatures, (Fig. 7). Upon cooling from 610°C, the Curie point was lowered by 
about 60°C. Ten hours again elapsed before the final measurement was made, and 
the changes occurring between 62 and 20°C were of the same order of magnitude as 
those observed for F1 in the corresponding temperature range. 

(c) To determine the time interval during which the large increases in magnet- 
ization occurred the experiment was repeated with a third specimen, F3. Again this 
was heated above 600°C. After cooling to 50°C, it was taken out of the furnace and 
transferred to an apparatus where volume susceptibility & was measured in a small 


field, at repeated intervals. 


Table Il. Specimen F3—Volume susceptibility 
t=time in minutes since first reading was taken. 
k=incremental volume susceptibility in c.g.s. units X10-8. 


t k t k t k 
—* | 4,700 7.5 19,170 15.0 21,270 
0 9,630. 8.5 | 19,680 16.0 | 21,460 
4.0 15,210 10.5 20,850 17.0 21,590 
6.0 17,600 13.5 | 21,170 19.0 21,600 


* Reading taken before the specimen was heated. 


Table II shows that k increased steadily for the first 10 minutes and then 
approached a stable value. The reading at t=19 minutes should correspond to a temper- 
ature of 20°C and at ¢=0 the temperature should be in the range 30-40°C. In the 
time elapsed & had increased by a factor of 2.2, which is of the same order as the 
increases in J,’ for Fl and F2, in roughly the same temperature interval. The con- 
clusion appears warranted that temperature rather than time effects were responsible 
for the peculiar magnetization changes observed in specimens F1-F3. It may be 
suggested that the first heating caused chemical alterations in the material, resulting 
in the formation of a new constituent having its Curie point near room temperature, 


aoe ‘y 


or slightly above. 
(d) To investigate the possibility of chemical action, a further specimen, /4, 


was subjected to the same heat treatment as F1-F3 and then cooled to —50°C by 
means of solid carbon dioxide. Inspection of Fig. 8, which shows the curves obtained 
in a field of 1200 Oe., combined with those previously obtained for F1, but adjusted 
to the new scale, reveals an intensely magnetic constituent with a definite Curie point 
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in the low temperature range. Independent estimates from the J/,’-T and J/,/-T curves 
gave 15°C and 22°C, respectively. This compares with the value of 61°C, quoted by 
Wologdine, [19] the only reference to the Curie point of franklinite discovered in the 
literature. The coercive force minimum near 20°C is perhaps due to the interaction 
of the two constituents. A similar experiment with a specimen from the original 
franklinite sample, but which had not been heated, failed to disclose a low-temperature 
constituent of major intensity. The phenomena observed can thus be satisfactorily 
explained on the basis of chemical action caused by the heat treatment. 

G.D. Nicholls, of the Department of Geology, University of Manchester, has 
examined some of the above specimens in powder form. He found that the original 
“franklinite ” sample was really a mixture of 

1. An opaque magnetic constituent, probably franklinite ; 

2. A green mineral, believed to be a pyroxene, or possibly willemite (zinc 
silicate) ; 

3. A reddish brown mineral, perhaps tephroite (manganese silicate). 

After heat treatment, mainly the opaque mineral, but also some of the silicate 
remained, thus confirming that chemical change had taken place. An X-ray powder 
photograph showed the cell dimensions of the opaque mineral to be very similar to 
those of magnetite, as would be expected in franklinite, but no detailed measurements 
were possible in this case. 


It may be plausibly suggested that the chemical changes caused by heating P 


resulted in a transfer of some Zn and Mn into the crystal lattice of the franklinite, 


which is a mixed ferrite, containing variable proportions of these atoms in its natural 


state already. It has been shown [20] that the addition to a ferromagnetic ferrite of 
zinc ferrite, which is paramagnetic at normal temperatures, produces mixed crystals 
_ whose Curie point is below that of the ferromagnetic ferrite. If manganese ferrite, 
whose Curie point is at 510°C, is also added, the magnetic properties of the ferro- 
magnetic ferrite, in this case franklinite, are probably further modified. A mechanism 
on these lines may account for the phenomena observed with specimens EE -F4, 


a oer 


dO} ThE, pegs atier: was extended to ‘specimen 5, cut from Cee new 
sample of frankl: 48 to 


eee 
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concerned. 


A more precise explanation of 


the phenomena here observed with 


° 200 400 600 ° 200 400 600 


franklinite would require a good deal 


of further work, including a detailed 


study of the chemical processes occur- 
© 200 400 600 °? 200 400 600 : i 
ring during heat treatment. 
(4) Basalt 
Basalts commonly contain 


magnetite or titanomagnetite as minor 


constituents, in which case they exhibit 


: 
© 200 400600 ° 200 400 600 200 400 600 


Fig. 8 Specimens Fig. 9 Specimen Fig. 10 Specimen ferromagnetic properties. In his work 


F4 and F1. FS. Franklinite N.‘‘ Normal’? on the dykes and lava flows of Mull, 


Franklinite S : ; ; 
: : Basalt. Vincenz [2], [3] measured the directions 
Points obtained © before heating ; @ whilst heating ; 
© whilst cooling. 
H, and H,! in oersteds; I,, I,’ I, and I,! in arbitrary @d samples from each rock formation 


units. he investigated, and computed a mean 


of magnetization of a number of orient- 


direction for that formation. The significance of this mean was assessed by a statistical 
analysis. 

The experiments described below have been conducted principally to compare 
the temperature curves of rocks whose natural thermoremanence had been respectively 
normal and inverse in direction. To be fully representative of these two types of 
behaviour, the specimens were chosen from rock formations whose computed mean 
direction of magnetization had been shown to be significant. The normally polarized 
rock was taken from a “normal” lava flow and the rock with inverse magnetization 
from a dyke whose polarization was inverse. To permit accurate measurements, the 
specimens were selected from a group having reasonably large susceptibilities when 
measured in weak fields (k>4000x 10-8 c.g.s. units). Finally, neither of the specimens 
used had undergone heat treatment prior to the present investigation. 

Both specimens were examined at H)=2300 Oe., as they cooled from 600°C. To 
facilitate comparison between them, J, and J, were plotted on a scale which made 
their values at 20°C equal to 1000. 

(a) The “normal” basalt, specimen N, came from a lava flow in the north eastern 
part of Mull. The average value of the Curie point obtained from the three curves 
(Fig. 10) was 572°C, or less than 10°C below that of pure magnetite. The pronounced 
bend in the H,.-T and J,-T curves above 250°C suggests the presence of a second 


_- 


constituent, with its Curie point near 300°C. z 
Koenigsberger [14] examined four German late Tertiary or Quaternary basalts, 


containing magnetite along with titanium dioxide and excess ferrous oxide and, in one 
case, pyrrhotite as well. All four specimens had two Curie points each, at temperatures 
ranging from 450 to 620°C, and 170 to 300°C, respectively. H, of the basalts at room 
temperature was between —90 and —175 Oe., compared with —354 Oe. obtained after 
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heat treatment in the present case. However, the coercivity of the high Curie point 
constituent in N was —188 Oe. at 350°C, where the H,-T curve is nearly horizontal ; 
this value appears to be more in agreement Koenigsberger’s results. The large 
coercivity (—534 Oe.) before heat treatment may be due to a very small grain size of 
the ferromagnetic constituents, crystals in volcanic rocks being commonly tiny. The 
changes in J, and J, produced by heat treatment were small, amounting to —10% — 
and +3%, respectively. 

(b) The “reversed” specimen, R, was a dolerite basalt from a dyke near 
Dishig, on the south shore of Loch na Keal, Mull. The main features of the three 
curves (Fig. 11) are simple, but it appears that a major Curie point at about 500°C 
is followed by that of a minor constituent, at about 600°C. Since /, and J, are very 
small above 500°C, the trend of H, is doubtful in this region, and it may vanish 
between 540 and 600°C. 

A major second constituent with its Curie point in a lower range of T would 
have been required, for example, by Néel’s two-component mechanisms of reversal. 
Apart from the small irregularities observed, however, specimen R shows the behaviour 
of a relatively pure ferromagnetic material, whilst specimen N, which would be 
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caused no important chemical changes and that conclusions drawn from the laboratory 
data are here applicable to the conditions under which the dyke cooled through its 
Curie point. 

J. Singh [21], until recently, of the Geophysics Department, Imperial College, 
investigated the loss of natural remanent magnetism on heating and examined an 
inversely magnetized specimen frém the same flow as specimen R. His results are 
shown in Fig. 11b’, which bears a striking resemblance to the -T curve of specimen 
R. Again, no evidence of a second constituent in significant proportions was found. 

It may be suggested that the basalt curves obtained here are, apart from the 
second constituent of specimen AN, more truly representative of the properties of 
magnetite or titanomagnetite than those of the magnetite specimens themselves, M1 
and M2 having changed completely after heating, whilst M3 was not examined under 
saturation conditions. 

From the geophysical point of view, the essential feature of both H,-T curves 
discussed in this section appears to be the fact that high H, values were attained in 
the range 400-500°C. This implies that the basalts were capable of retaining thermo- 


remanence, even at high temperatures, in opposition to an acting field. Together with 


the failure to find evidence in favour of the two-component mechanism, the data 
therefore tend to support the hypothesis that a past reversal of the earth’s magnetic 
field has been responsible for the inverse magnetization of the rock formations con- 
cerned. However, since they were based on two specimens only, the results must be 
regarded as supplementary to other evidence. 


Summary 


Most of the minerals and rocks examined became fully or nearly saturated in 
an external field of 23000e. The specimens behaved like normal ferromagnetic 
materials, their 7,-T curves broadly resembling those of the ferromagnetic elements. 
H, was always large at normal temperatures and, in most cases, it still exceeded 
—100 Oe. at 500°C. It is safe to conclude that the rocks and minerals concerned are 
capable of retaining thermoremanence for long periods, even at high temperatures. 

Irregular features in the temperature curves of some specimens indicated the 
presence of more than one ferromagnetic constituent. Probably due to chemical 
action in the material, the hysteresis loops measured at 20°C before and after the heat 
treatment differed considerably in some instances. In the case of a certain franklinite 
with its Curie point at about 600°C a new, strongly magnetic, constituent with its 


Curie point near room temperature appeared after heating. —— 
A comparison of data from two basalts, with respectively normal and inverse 


- directions of natural thermoremanence, gave no indication that a property of the 


material caused the inverse polarization in one of the rocks. The results supported 
the hypothesis that such magnetizations owe their direction to past reversals of the 


geomagnetic field, a 
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